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ABSTRACT 

It is expected that energetic solar flares releasing large amount of energy at 
the photosphere may be able to excite the acoustic (p-) modes of oscillations. We 
have determined the characteristic properties of mode parameters by applying 
the ring diagram technique to 3-D power spectra obtained for solar active region 
NOAA 10486 during the long duration energetic X17.2/4B flare of October 28, 
2003. Strong evidence of substantial increase in mode amplitude and systematic 
variations in sub-surface flows, i.e., meridional and zonal components of velocity, 
kinetic helicity, vorticity, is found from comparison of the pre- to the post-flare 
phases. 

Subject headings: Sun: activity — Sun: flares — Sun: helioseismology 



1. Introduction 



Leighton et al.l (119621 ) discovere d solar 5-min ute oscillations while measur ing Doppler 



shifts in photospheric spectral lines. lUlrichl (jl970l ) and iLeibacher fc Steiru (119711 ) later inter- 
preted these oscillations to be due to acoustic (or p-) modes trapped in the sub-photospheric 
layers. Accurate determination of these modes provides a powerful diagnostic tool for prob- 
ing the solar interior. It is generally believed that these modes may be either intrinsically 
overstable or stochastically excited by turbulent convection. 



Characteristics of the surface p-modes and sub-surface flows are essentially described by 
shorter wavelength modes that are trapped near the surface. These modes can be studied 
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using the ring diagram analysis which uses three dimensional data cubes (16° x 16° x 1664), 
first two correspond to spatial size of the active region and the third is time in minutes (IHill 
19881 ). Other techniques of local helioseismology, generally used for studying various as- 
pects of sub-photospheric characteristics of active re gion interiors, are t ime-distance analysis 
(IZhao fc Kosovichevll2004j ). and seismic holography (IBraun et al.l 120041 ) . 



Wolflg (I1972I ) first suggested that energetic flares may be able to excite acoustic waves by 



exerting mechanical impulse of the thermal expansion of the flare on the photosphere. They 
estimated the damping times to be longer than a day for the free modes. Subsequently, active 



region (AR) related effects on p-mode parameters have been reported 
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Donea et al.l (119991 ) have also reported excitation of fiare-related waves on the solar surface, 
however, these pertain to traveling waves as opposed to standing waves which constitute the 
normal modes of solar oscillations. 

The spatial extent of flares is usually much smaller than the spatial sizes of 16° x 16° used 
in GONG and SOHO/MDI data-cubes for ring diagram analysis. Also, any mode amplifica- 
tion induced by transient flare events has to essentially compete with the absorption effects 
associated with intense magnetic fields of sunspots. Therefore it is expected to be rather 
difficult to conclusively resolve flare-related effects by averaging techniques. Here one can 
estimate the damping time scale of p-modes from the observed mode width, generally in the 
range of 15-100 ^Hz for high-degree modes. This gives a lifetime in the range of 0.5-3 hr 
assuming that the width is mainly due to the finite life of the modes. Acoustic waves travel 
outwards from the flare site with the speed of sound varying in the range from 10 to 50 
km s"^ for the depths in which the high-degree modes are trapped. Thus, over the damping 
time of the order of an hour, these waves would travel a distance from 36 to 180 Mm, i.e., 
comparable to the spatial extent of the active region. A large fraction of energy dissipa- 
tion occurs within this region if the flare is temporally located well within the data-cube. 
Therefore, it is expected that flare-related effects would last over several hours and should 
be detectable in the temporal averages carried out over a day as in the ring analysis. On the 
other hand, if the flare occurred around the end time of the data-cube, the flare effects would 
not be detectable as its effects would be contained mostly in the subsequent data-cube. 

Relation between photospheric motions and flare activity has been studied earlier. For 
example, correlation tracking studies have discovered a few cases of sma ll-scale vorticity a t 
the surface preceding a flaring event with timescales around an hour (lYang et al.l 120041). 



Some cases of sunspot rotation in flaring active regions have also been found (IBrown et al. 



2OO3I ). Relationship of flare activity with some statistical properties of active regions has 
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been studied bylLu| (Il995l ) and lAbramenko &: Longcopd (120051 ) . among others. It is believed 
that twisted flows in the sub-surface layers cause foot-point motions of sunspots observed at 
the photospheric layer leading to unstable magnetic topologies in the overlying regions. 

Although the average properties of p-mode parameter s and sub-surface flows have been 
studied with magnetic /flar e activities by many r esearchers ( Ambastha et al.ll2004 : Howe et al. 
20041 : iKomm et allboOSbl lal: | Mason et al.ll2006l ). the change of p-mode parameters and sub- 
surface flows during flares is not well understood. Some efforts have been made to examine 
the difference betwe ei i flare productive and quiet re gions to infer any activity related effects 
(IMason et al.ll2006l ). lAmbastha et al.l (l2003l . 120041 ) have found that power in p-modes ap- 
pears to be larger during the period of high flare activity as compared to that in non-flaring 
regions of similar magnetic field strength. However, the pre-, peak and post-flare signatures 
of flares on jo-mode parameters and sub-surface flows still require careful analysis. We have 
studied the long duration, energetic X17.2/4B flare of 28 October 2003 in NOAA 10486. We 
found strong evidence of flare-related changes in p-mode parameters and sub-surface flows 
by placing the flare at different temporal positions in appropriately constructed data-cubes. 



2. The White Light Flare of NOAA 10486 

In order to determine any flare-related effects on p-mode parameters, we have considered 
one of the most energetic events that occurred in NOAA 10486, i.e., the white light super- 
flare of 28 October 2003, classified as X17.2/4B event. During the period of the data-sets 
used in our study, some other C and M class flares also occurred in NOAA 10486, but 
the X17.2/4B super flare was the dominant event during the 24-hr period. According to 
Geostationary Operations Environmental Satellite (GOES) X-ray observations, this flare 
started at 09:51 UT, reached the maximum phase at 11:10 UT, and decayed at 11:24 UT, i.e., 
lasted over more than 90 minutes (Figure [1]). However, it is evident that the integrated X-ray 
flux remained at a very high level well beyond 11:24 UT. Even if a background corresponding 
to the Ml level (10^^ watts m~^) is considered, the X-ray flux gradually reduced to this 
level only around 16:00 UT. In Ha also, it is reported to have lasted for more than four 
hours. Therefore, this long duration event (LDE) was a particularly well suited case for an 
investigation of flare-related helioseismic effects. 
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3. The Data and Analysis 



We have used high resolution GONG Dopplergrams obtained at one minute cadence. 
The data-cubes have 1664 minutes' duration and cover 16°xl6° area centered at 285° Car- 
rington longitude and -22.5° latitude covering NOAA 10486. The choice of area makes a 
compromise between spatial resolution, range of depths and resolution in spatial wavenum- 
ber in the power spectra. A larger size will allow access to the deeper sub-surface layers, but 
only with coarser spatial resolution. On the other hand, a smaller size will not allow access 
to the deeper layers and also render the fitting of rings more difficult. The AR was located 
close to the disc center at the time of this flare, therefore, projection effects did not pose any 
serious difficulty in the analysis. 

The 3-D Fourier transform of the data-cube gives a trumpet like structure in the fre- 
quency domain [kx, ky, uj). A slice at fixed frequency u of this structure gives power concen- 
trated in concentric rings (jHilllll988l ). These rings provide the characteristic properties of 
p-modes and subsurface flows. The phase velocity of the acoustic wave is augmented by flow 
velocity (U) causing change in frequencies {Au = k.U) that perturbs the center and shape 
of the rings. Surface flows are derived from this perturbation using ring fltting. The p-modes 
of different wavelengths are trapped at different depths beneath the surface. Therefore, flows 
derived at the surface are weighted averages over depths. Inverting these characteristics of 
the modes give flows in the interior. 

As illustrated in Figure [H we constructed flve data-sets of 1664 minutes duration each 
with different starting time, such that the flare onset was placed in the beginning (-R5), one- 
fourth (-R4), center (-R3), three-fourth (_R2) and end (-Ri) of the data-cube's time-line. The 
16° patch consists of 128 pixels, giving a spatial resolution Aa; = 1.5184 Mm, i.e., the k- 
number resolution, Ak = 3.2328 x 10^^ Mm~^, and a Nyquist value for the harmonic degree, 
/ = 1440. The corresponding range in (^kx, ky^ space is from -2.069 Mm"^ to 2.069 Mm'^ 
The temporal cadence and duration of data-sets give Nyquist frequency of 8333 /iHz and 
frequency resolution of 10 //Hz, respectively. 

To determine the surface mode par ameters of jp-modes , we have carried out ring- diagram 

analysis using the den se-pack technique (IHaber et al.ll2002l ) adapted for GONG data fICorbard et al, 



2003 



Hill et al.ll2003l ). Images were remapped around the central position (285°, -22.5°) using 



transverse cylindrical projection to obtain equi-distance spatial sampling interval required 
for Fourier transform. To remove the effect of differen tial rotation, the remapped images 
were tracked with differential rotation rate of the Sun (jSnodgrassI Il984j ) . The tracked im- 
age cubes are apodized before being Fourier transformed. The main ring-fltted parameters 
(i.e., surface acoustic mode parameters) determined in our study are radial order (n), degree 
(/), mode amplitude (A), mode width (F), zonal velocity (Ux) and meridional velocity (Uy). 
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The commo n mode parameters in all sets were corrected for filling factor using a method 



described in lKomm et al.l (120001 ) . Zonal {u^) and meridional (uy) components of sub-surface 



flows were d erived by regularized least square (RLS) inversion from the surface to a depth 
of ~20 Mm Jxhompson et aljll996l : Iflaber et al.1l2002h . 



The topology of fluid is measured by kinetic helicity. It is defined a s the volume integral 
of th e dot product of the velocity (u) and vorticity (u) of the flows (IMoffatt fc Tsinober 
1992h : 



Hi< 



u ■ uj dV 



V 



where, a; = V x u is the vorticity vector. It is a measure of circulations per unit area, i.e., 
twist of the flow, while its dot product with velocity gives variation of twist with flows. We 
estimated the vertical compo nent of velocity u^, fr om the divergence of horizontal components 
assuming mass conservation (IKomm et al.ll2005al ). Now, having all three components of flows 
with depths, we can derive kinetic helicity. As we do not have the flows at each point of the 
active region from ring analysis, we can derive only the average helicity over the entire area, 
termed as helicity density {hx = u • V x u). The numerical derivatives of flows are derived 
using Lagrange three point interpolation. 



4. Results and discussions 

Mode parameters obtained for the five data-sets are shown in Figure [T]for comparison of 
the flare-related effects. Surface acoustic mode parameters and the corresponding inverted 
parameters for these data-sets are shown in Figures [2H11 

Figure [2] shows that relative difference in mode amphtude obtained for i?2 to R5 where Ri 
is used as the reference. This quantity is found to be the largest for as compared to that for 
other sets. The difference between the amplitudes increased with frequency i/ > 2000 /iHz, 
and reduced with increasing radial order from n = to 5. The large difference of mode 
amplitude at lower radial orders may be attributed to the flare effect as these modes are 
confined near the surface and hence, more likely to be affected by the surface activity. 

The flare was placed at the beginning in R^, hence the overall helioseismic contribution 
of this LDE event, extending from the pre- to post-flare phase, is expected to remain entirely 
within this data-cube. Correspondingly, the mode amplitude was found to be the largest in 
i?5 as expected. On the other hand, it was smallest in Ri where the flare was placed at the 
end of the data-cube, therefore post-flare effects were not covered by this data-set. From the 
systematic trend of variations in the mode amplitude from R^ to Ri, it is evident that the 
energetic flare indeed gave rise to a significant amplification of mode amplitude. 
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Mode width is also found to be large at all frequencies and radial orders for R^. The 
changes ar e found to be l arger at lo wer frequencies i n aplyin g that higher-n modes have shorter 
life times. iDuvall et al.l (119881 ) and iBurtseva et al.l (120091 ) have also found that mode width 
increased with frequency as well as with degree of the modes. However, care is needed in 
this interpretation as significant contribution to mode width may arise due to the available 
limited resolutions in wavenumber and frequency. 

Using the ring diagram fits, we have calculated the horizontal velocities in order to 
check if there is any variation in surface velocities due to flares. We found large variations 
in surface zonal and meridional velocities at all frequencies for the radial orders n = to 
5 from i?i to R5. Also, the deviation for both the components of velocity is larger for 
as compared to other data-sets indicating the variation in flows from the pre- to post-fiare 
phases. 

The fitted velocities for each mode can be inverted to calculate the zonal (ux) and 
meridional (uy) components of velocity as a function of depth. These were derived by regu 



larized least square (RL S) inversion from the surface to a depth of 20 Mm (IThompson et al. 



19961 : iHaber et al.ll2002l ). Both and Uy exhibited significant systematic changes with the 
flare. The profiles of zonal velocities Ux with depth are shown for Ri to R5 in Figure [3] {left 
panel) . A large decrease in zonal velocity around the depth ^ 4 Mm is common in all cases. 
However the depth where the minimum of Ux occurred varied from the shallowest for Ri 
to the deepest for R^. Figure [3] {right panel) shows the coresponding changes in meridional 
velocities Uy for Ri to R^. It is found that Uy decreased rapidly from the surface to the depth 
of ~ 1 Mm for all the data-sets. Thereafter, it increased and attained a peak at ~ 6 Mm. 
It is inferred that the gradient in Uy is largest for Ri (i.e., before the flare) and it decreased 
systematically from R2 to R^ as the post-fla re effects increased. Th is confirms the earlier 
reports about meridional velocity gradient by lAmbastha et al.l (120041 ). 



Figure H] shows the profiles with depth of (a) divergence of horizontal components of 
velocity, (b) vertical component of velocity Uz, (c) vertical vorticity and (d) helicity density. 
Interestingly, it is observed that flare-related variations in helicity density hk are confined 
within a shallow region from the surface to a depth of around 3 Mm. For depths greater 
than 3 Mm, hk derived for all the five cases approached toward the same value. 

Changes in the topology of sub-surface flows during pre- to post-flare phases are illus- 
trated in Figure H^-d. In Ri, the divergence of horizontal component {uh) near a depth 
of ~ 1 Mm changed from a small negative to a large positive value for R2, and decreased 
from to -R5 (Figure S^). This suggests that the perturbation in the flows at this depth 
gradually returned back to the initial or pre-flare state after the flare decayed. On the other 
hand, the large positive divergence observed at depths below 3 Mm decreaseed in magnitude 
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from Ri to R^. 

The vertical component of flows (uz) changed from upward direction in Ri to downward 
direction in R2 at a depth of ~ 1 Mm, and then back to the upward direction in R^ through 
-R3 - -R4 (Figure Hb). The Uz profiles converge to the same negative value (~ —0.1 ms~^) at 
the depth of ~3 Mm, and then diverges systematically from Ri to R5 as the depth increases. 
The vertical component of vorticity vector (uz) is shown in Figure Ht for Ri to -R5. The 
peak in Uz near 2 Mm and the trough near 6 Mm manifest the bipolar structure of the AR. 
Kinetic helicity (hx) changed from a small positive value for Ri to a large negative value 
for R2, increasing toward the level of Ri for R^ - R^ (Figure HJi). However, it remained 
nearly constant for all the data-sets at depths below ~3 Mm. These systematic changes in 
sub-surface parameters {uh, Uz, oJz and hx) from Ri to R^ provide unambiguous evidence of 
the relationship of the large flare with sub-surface dynamics of the active region. 



5. Conclusions 

This study of an extremely energetic and long duration X17.2/4B flare in NOAA 
10486 gives a clear indication of fiare-related variations in the p-mode parameters and sub- 
photospheric flows in NOAA 10486. Changing the temporal position of the flare within the 
Doppler data-cubes obtained for the active region amounts to the changing level of the pre- 
and post- effects in the data-set. The ring diagram analysis of the different data-sets thus 
constructed provides the following important results: 

(i) The amplitude of p-modes increased up to 150% in the case of the flare placed in 
the beginning of the data-set as compared to the case when flare was placed near the end or 
outside the data-cube. A similar result is obtained for p-mode energy also as expected due 
to its relation to the amplitude, manifesting the rate of energy supplied to the p-modes by 
the flare. 

(ii) Amplitude and energy of the modes decreased with radial order indicating that the 
effect of the flare decreased with increasing depth. Furthermore, we found that the amplitude 
and energy of modes increased with frequency. This suggests that modes with high natural 
frequencies are amplified more by the flare as compared to the low frequency modes. 

(iii) The gradient in meridional velocity observed at depths 2-6 Mm decreased with the 
flare which manifests relationship of the flare activity with the sub-surface flows. 

(iv) The vertical component of the flow changed to the downward direction during the 
flare and then returned back to the pre-flare state. This is an evidence of downward moving 



-8- 



material during the flare. 

(v) Divergence of the horizontal component of flow near the surface changed from neg- 
ative to positive and then back to the pre-flare state after the flare. 

(vi) The sub-surface flow possessed a bipolar structure with one pole located at the 
depth of around 2 Mm and the other at 6 Mm. Near the surface, the flow was twisted during 
the pre-flare phase, which relaxed after the flare. 

(vii) The kinetic helicity of sub-surface flow around the depth of 1 Mm changed from 
positive to negative during the pre- to the peak phases, and returned back to the pre-flare 
level after the post-flare phase. 

In summary, this study provides strong evidence about the role of a large flare in modi- 
fying various acoustic mode parameters. Depending upon the temporal location of the large 
flare in NOAA 10486, systematic variations in mode parameters have been found for dif- 
ferent data-sets with varying levels of pre- and post-flare effects. These results also suggest 
that flare-related changes in the acoustic mode parameters are detectable by ring diagram 
analysis provided that the pre- and post-flare phases of the flare (or flares) are well covered 
by the data-cube. Finally, we conclude that the sub-surface flow topology can be used as a 
proxy for flare-forecasting. 

This work utilizes data obtained by the GONG program operated by AURA, Inc. and 

managed by the National Solar Observatory under a cooperative agreement with the National 
Science Foundation, U.S.A. The integrated X-ray flux data was obtained from GOES which is 
operated by National Oceanic and Atmospheric Administration, U.S.A. The authors would 
like to thank H.M. Antia for many useful discussions during the course of this work and 
providing mode inertia table for computing mode energy. 
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Fig. 1.— GOES-12 integrated X-ray flux for 1.0 - 
line). The horizontal lines with labels Ri,i — 1, . . . 
for constructing the flve data-sets. 



O.SA (sohd line) and 0.5 - 0.4A (dotted 
5, represent different time periods taken 
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Fig. 2. — Relative difference in mode amplitude for the data-sets Ri,i — 2, . . .5, for radial 
orders, n — - 5. Here, Ri is used as the reference. 
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Fig. 3. — Zonal (left panel) and meridional Uy (right panel) velocity profiles with depth 
from to 12 Mm for the data-sets Ri, i = 1, ... 5. 
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Fig. 4. — Inverted profiles of (a) divergence of horizontal components of velocity, (b) vertical 
component of velocity Uz, (c) vorticity and (d) helicity density, with depth from to 12 Mm 
for the data-sets Ri,i — 1, . . .5. 



